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Abstract 

Deep-inelastic e^p scattering at high squared momentum transfer Q 2 up to 30000 GeV 2 is 
used to search for eq contact interactions associated to scales far beyond the Hera centre 
of mass energy. The neutral current cross section measurements da/dQ 2 , corresponding 
to integrated luminosities of 16.4 pb" 1 of e~p data and 100.8 pb _1 of e + p data, are well 
described by the Standard Model and are analysed to set constraints on new phenomena. 
For conventional contact interactions lower limits are set on compositeness scales A ranging 
between 1.6 — 5.5 TeV. Couplings and masses of leptoquarks and squarks in imparity 
violating supersymmetry are constrained to M/X > 0.3 — 1.4 TeV. A search for low 
scale quantum gravity effects in models with large extra dimensions provides limits on the 
effective Planck scale of Ms > 0.8 TeV. A form factor analysis yields a bound on the 
radius of light quarks of R q < 1.0 • 10~ 18 m. 
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1 Introduction 



Deep inelastic neutral current scattering ep — > eX at very high squared momentum transfer Q 2 
allows one to study the structure of eq interactions at short distances and to search for new 
phenomena beyond the Standard Model (SM). The concept of four-fermion contact interactions 
(CI) provides a convenient method to investigate the interference of any new particle field as- 
sociated to large scales with the 7 and Z fields of the Standard Model. This paper considers 
conventional contact interactions, such as general models of compositeness and the exchange 
of heavy leptoquarks and supersymmetric quarks, as well as low scale quantum gravity effects, 
which may be mediated via gravitons coupling to Standard Model particles and propagating 
into large extra spatial dimensions. 

The present analysis is a continuation of previous studies [1] based on e+p data [2]. New 
HI cross section data on e~p scattering [3] and on e + p scattering data [4] lead to a substantially 
higher sensitivity to new physics phenomena. Similar studies of eq contact interactions have 
been performed by other experiments at Hera [5], Lep [6] and Tevatron [7], providing 
results comparable to those presented here. 



2 Data and analysis method 

The data have been collected with the HI detector at Hera and correspond altogether to an 
integrated luminosity of C int = 117.2 pb _1 . They consist of three data sets of different lepton 
charge and centre of mass energy y/s: 

reaction C int [pb _1 ] y/s [GeV] run period ref. 



e + p -> e+X 35.6 ± 0.53 301 1994 - 1997 [2] 
e'p^e-X 16.4 ±0.3 319 1998 - 1999 [3] 
e+p -> e+X 65.2 ±0.95 319 1999 -2000 [4] 

The cross section measurements da/dQ 2 extend over a large range of squared momentum trans- 
fers 200 GeV 2 < Q 2 < 30000 GeV 2 for inelasticity y < 0.9. Details can be found in the quoted 
references. In this paper a contact interaction analysis is presented for each of the two recent 
data sets and the combined data including the previous measurements. 

The phenomenological models under study and the analysis method are described in more 
detail in ref. [1]. The analysis investigates the measured cross sections da/dQ 2 and performs 
quantitative tests of the SM or CI models, applying a minimisation of the x 2 function 



Here a t cxp and af 1 are the experimental and theoretical cross sections for the measurement point 
% and A<7j is the corresponding error including statistical and uncorrelated systematic errors 
added in quadrature. The functions A ik (£ k ) describe correlated systematic errors for point % 
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associated to a source k. They depend on the fit parameters e k , which are effectively pulls 
caused by systematics. In general the influence of the correlated systematic errors is small. The 
following dominant sources of correlations are included: an overall normalisation uncertainty 
of 1.5 — 1.8 % depending on the run period, experimental uncertainties on the scattered lepton 
energy and angle and an uncertainty on the strong coupling a s (M z ) = 0.118 ± 0.003 entering 
the SM prediction. In a combined x 2 analysis using eq. (1) all data sets are treated as indepen- 
dent samples with individual normalisation and measurement uncertainties. The influence on 
the fit results of correlations between data sets is negligible. 

The cross sections da(e~p — > e~X)/dQ 2 [3] and da(e + p — > e + X)/dQ 2 [4] from the new 
data are shown in figure l. 1 They are very well described by the Standard Model expectations 
over the full Q 2 range, over which they vary by six orders of magnitude. Cross section calcu- 
lations in the Standard Model are performed in the DIS scheme in next-to-leading order (NLO) 
QCD using the corresponding CTEQ5D parton parametrisation [8]. These parton distributions 
were obtained including the small amount of 2.7 pb _1 of 1994 HI data with Q 2 < 5000 GeV 2 
only and are thus essentially uncorrelated with the present data. 

Fits to the SM prediction yield x 2 — 10-4 for 17 degrees of freedom (dof) with a fitted 
normalisation of 0.998 for the e~p data and x 2 /dof = 13.9/17 with a normalisation of 1.006 
for the e + p data, using statistical errors from the experiment and including all systematics as in 
eq. (1). Applying other parton density functions like CTEQ6 [9], MRST [10] or GRV [11] leads 
to slight changes of the normalisation within errors, but yields equally good agreement with the 
data. 

The cross section measurements do not show significant deviations from the Standard Model 
and will be used to set limits on couplings from processes mediated through contact interactions. 
Since the concept of contact interactions is an effective theory, theoretical expectations cannot 
be formulated consistently in NLO. A sensible approach is to reweight at each Q 2 value the SM 
NLO cross section by of°(SM + CI)/<3f °(SM), the ratio of leading order (LO) cross sections 
with and without inclusion of contact interactions. The x 2 function of eq. (1) can be applied to 
evaluate the sensitivity of the data to a certain CI scenario and to determine its parameters. The 
corresponding limits at 95% confidence level (CL) are derived by using a frequentist approach 
as described in the Appendix. Systematic uncertainties are included in this procedure. This 
determination of limits is different from the method used in our previous publication [1]. 

3 Contact interaction formalism 

The most general chiral invariant Lagrangian for neutral current vector-like four-fermion con- 
tact interactions can be written in the form [12, 13] 

£v = V 9 a b (e a l^e a )(q b Yq b ) . (2) 

q a,b = L,R 

For lepton e and each up-type and down-type quark flavour q with the corresponding currents e a 
and q b there are four coupling coefficients r] q ab = e ab (fl'/A^) 2 , where a and b indicate the L (left- 
handed) and R (right-handed) fermion helicities, g is the overall coupling strength, A q ah is a scale 

'Corresponding figures for the e + p cross sections at = 301 GeV [2] are shown in [1]. 
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parameter and e ab determines the interference sign with respect to the Standard Model currents. 
Any particular model, such as compositeness or the exchange of leptoquarks or supersymmetric 
quarks, can be constructed by an appropriate choice of the couplings r] q ab . The phenomenological 
models of interest and their analytical treatment are discussed in more detail in [1]. 



4 Compositeness scales 

In general models allowing for fermion compositeness or substructure it is convenient to choose 
a coupling strength of g 2 = An and to assume a universal scale A for all quarks. The contact 
interaction coefficients then read 

q _ 47T 
Vab ~ € ab ■ 

Various scenarios of chiral structures, e.g. pure L(eft) and R(ight) or V(ector) and A(xial vector) 
couplings, are defined by setting particular chiral contributions to values of e a b = ±1 and setting 
all other combinations to zero. 

The sensitivity to CI models is tested by determining the quantities e/A 2 in a x 2 fit using 
the experimental statistical and systematic errors and leaving the sign of interference free. The 
results from a combined fit to all data sets are shown in figure 2. The data tend to prefer negative 
interference, but are consistent with the Standard Model within two standard deviations in all 
scenarios 

Lower limits on compositeness scale parameters A ± , associated to positive or negative in- 
terference, are derived from a frequentist approach and are summarised in table 1. Despite 
the lower integrated luminosity, the e~p data exhibit for some models comparable (LL and 
RR couplings) or even higher (AA + coupling) sensitivity than the e + p data. The two lepton 
charges complement each other and a combined analysis of all e ± p data sets yields substan- 
tial improvements in most scenarios compared with the previous publication [1]. The results 
are presented in figure 3. The lower limits on compositeness scale parameters vary between 
1.6 TeV and 5.5 TeV depending on the chiral structure. Choosing different parton distribu- 
tions [9-11] changes the quoted limits typically by a few per cent, at most by 15 %. The most 
restrictive bounds are observed for the VV model, where all chiral components contribute with 
the same sign. As an illustration, examples of fits to the e~p and e + p cross sections are shown 
in figure 4. 



5 Leptoquarks 

Leptoquarks couple to lepton-quark pairs and appear in extensions of the Standard Model which 
try to connect the lepton and quark sectors. They are colour triplet scalar or vector bosons, 
carrying lepton (L) and baryon (B) number and a fermion number F = L + 3 B. Since F = 2 
for e~q and F = for e + q states, one expects different sensitivities to particular leptoquark 
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types from electron and positron scattering. For high enough mass scales the leptoquark mass 
M L q and its coupling A are related to the contact interaction coefficients via 



'U tab M 2 q ■ 

Within the model of [14] the production and decay modes are fixed and the relative contributions 
e q ab have been calculated [15]. The notation, quantum numbers and couplings of the various 
leptoquarks are given in table 2. 

The analyses of the cross section measurements do not show an indication of a leptoquark 
signal. The results of fits for each type of leptoquark are interpreted in terms of limits on the 
ratio M L q/\ and are summarised in table 2. In some cases, e.g. and V X L , the e~p data 
give more restrictive bounds than e + p scattering despite the lower integrated luminosity. This 
sensitivity is illustrated in figure 4, which shows possible contributions of the leptoquarks S'f 
and to the e~p and e + p cross sections. The two leptoquarks differ by their spin and couple 
with the same chiral structure but different strength and sign to u and d quarks. This emphasises 
the complementary role and importance of both electron and positron beams. The combined 
analyses of all e ± p data further constrain the search for leptoquarks, reaching exclusion values 
of up to M L q /A = 1.4 TeV. Note that upper bounds on the coupling strength A can only be set 
for leptoquark masses exceeding the accessible centre of mass energy of Hera. 



6 Squarks in R p violating supersymmetry 

In the most general formulation of supersymmetry there exist operators which couple a lepton- 
quark pair to a squark, the scalar superpartner of a quark. Such couplings violate i?-parity (via 
lepton number violation), defined as R p = (_i) 3B + L + 2S ' with S being the spin. Thus R p = +1 
for SM particles and R p — — 1 for superpartners. This interaction allows single squarks to be 
produced or exchanged in deep inelastic scattering [16] via 

e + d R -> u L , c L , t L coupling A' yi , (3) 
e + u L -> d R , s R , b R coupling A' nfc . (4) 

The subscripts ijk of the coupling \' ijk describe the generation indices of the left-handed lep- 
tons, the left-handed quarks and the right-handed down-type quarks of the superfields, respec- 
tively. 

The e + q coupling of the left-handed up-type squark of reaction (3) is the same as that of the 
scalar leptoquark Sy 2 , and the coupling of the right-handed down-type squark of reaction (4) 
is the same as that of the scalar leptoquark Sq. Therefore the formalism and results of the 
leptoquark analysis can be directly applied. Limits on the ratio Mg/X' for R p violating squarks, 
assuming branching ratios Bq^ eq = 1, are given in table 3. Note that the squark generations 
cannot be distinguished in this analysis. 
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7 Large extra dimensions 



It has been suggested that the gravitational scale M s in 4 + n dimensional string theory may 
be as low as the electroweak scale of order TeV [17]. The relation to the Planck scale M P ~ 
10 19 GeV and the size R of the n compactified extra dimensions is given by Mp ~ R n M| +n . 
In some models with large extra dimensions the SM particles reside on a four-dimensional 
brane, while the spin 2 graviton propagates into the extra spatial dimensions and appears in 
the four-dimensional world as a tower of massive Kaluza- Klein states with a level spacing 
Am — 1/R. The gravitons couple to the SM particles via the energy-momentum tensor with 
a tiny strength given by the inverse Planck scale. However, the summation over the enormous 
number of Kaluza- Klein states up to the ultraviolet cut-off scale, taken as Ms, leads to an 
effective contact-type interaction [18] with a coupling 



The coefficient A depends on details of the theory and is expected to be of order unity. However, 
by convention, one also allows for a negative coupling and thus sets A = ±1. The scale de- 
pendence of gravitational effects is quite different from that of conventional contact interactions 
discussed in the previous sections. 

The cross section formulae for virtual graviton exchange in DIS have been calculated in [1] 
using the conventions and formalism of [18]. The interference of the graviton with the photon 
and Z fields has opposite sign for electron and positron scattering, as illustrated in figure 4. 
Lower limits on the scale parameter Ms, derived from fits to the da/dQ 2 distributions includ- 
ing gravitational contributions, are summarised in table 4. There is similar sensitivity to the 
effective gravitational scale for positive and negative interference, resulting in lower limits of 
M s > 0.82 TeV for A = +1 and M s > 0.78 TeV for A = -1. 

In other scenarios the gauge bosons 7 and Z are also assumed to propagate into extra di- 
mensions. This leads to analogous Kaluza-Klein states which couple to matter with electroweak 
strength and interfere with the ordinary gauge boson fields. In a specific model of 4 + 1 = 5 di- 
mensions [19] with compactification radius R = 1/M C , the sums over the exchange of Kaluza- 
Klein gauge bosons essentially modify the photon propagator 1/Q 2 — > l/Q 2 + tt 2 /(3Mq) 
and the Z propagator 1/(Q 2 + Mf ) -> 1/(Q 2 + Mf ) + tt 2 /(3 M 2 c ). A fit to the data yields 
M c > 1.0 TeV as a lower limit at 95% CL for the compactification scale. 



A fermion substructure can also be formulated by assigning a finite size to the electroweak 
charge distributions. It is convenient to introduce electron and quark form factors f(Q 2 ) which 
reduce the Standard Model cross section at high momentum transfer [20]: 



A 



71(3 ~ W s ■ 



8 Form factors 



f(Q 2 ) 



1 - 




(5) 
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da 



da SM 



feiQ 2 ) f q {Q 2 ) ■ 



dQ 2 



dQ 2 



(6) 
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Fits to the data yield upper limits on the particle size R = sj (r 2 ), taken as the root of the 
mean squared radius of the electroweak charge distribution. Assuming a point-like electron, i.e. 
setting f e = 1, the radius of the light u and d quarks can be constrained to R q < 1.0 ■ 1CT 18 m 
at 95% CL. If both electrons and quarks are assumed to have common form factors one obtains 
limits on fermion sizes of R eyQ < 0.7 • 10~ 18 m. The results of the analysis are given in table 5. 

9 Summary 

Neutral current deep inelastic e~p and e + p scattering cross section measurements are analysed 
to search for new phenomena mediated through (ee) (qq) contact interactions. The data are 
well described by the Standard Model expectations. The use of electrons and positrons pro- 
vides complementary information and a combined analysis based on an integrated luminosity 
of 117.2 pb _1 yields improved limits on scales of new physics. The present analysis supersedes 
previous results [1]. 

Lower limits at 95% CL on eq compositeness scale parameters A" 1 " are derived within a 
model independent analysis. They range between 1.6 TeV and 5.5 TeV depending on the 
chiral structure. A study of virtual leptoquark exchange yields lower limits on the ratio Mlq / A 
between 0.3 TeV and 1.4 TeV. Squarks in i?-parity violating supersymmetry with masses 
satisfying M^/X'^ < 0.43 TeV and Mj/A' nfc < 0.71 TeV can be excluded. Possible effects of 
low scale quantum gravity with gravitons propagating into extra spatial dimensions are searched 
for. Lower limits on the effective Planck scale Ms of 0.78 — 0.82 TeV are found. Allowing 
for Kaluza-Klein states of the SM gauge bosons results in a lower limit on the extra dimension 
compactification scale M c > 1.0 TeV. A form factor approach yields an upper limit on the size 
of light u and d quarks of R q < 1.0 • 10~ 18 m assuming point-like electrons. 
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Appendix Setting limits 

For the present analysis several methods have been studied to calculate limits and confidence 
levels and the final results quoted are based on the frequentist approach (see for example [21]). 
The x 2 function of eq. (1) is used as the quality measure of agreement between data and contact 
interaction models. It allows an easy implementation of systematic uncertainties. The other 
methods investigated differ in the definitions of the statistical error entering the total error Act; 
in the x 2 expression. 

Taking the statistical errors of the experiment, one can test the compatibility of the data 
with a certain model hypothesis and determine the corresponding parameter, for example 1/A 2 . 
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Each data point contributes with a fixed weight AcK~ 2 to the x 2 function, independent of the 
model parameter. However, downward fluctuations with respect to the model predictions get a 
larger weight (smaller statistical error) than upward fluctuations (larger statistical error). This 
property may enforce asymmetric situations and errors, as is observed in the present data, see 
figure 2. If the fitted parameter is found to be not significantly different from zero, the result 
can be converted into limits at a given confidence level (CL). A convenient way is to take the 
values corresponding to a change of x 2 (l/A ± 2 ) — Xsm = ^Xcl> f° r instance A%g 5% = 3.84 
for a 95% CL limit. This method was used in the previous publication [1] and is illustrated for 
several compositeness models by the x 2 distributions in figure 5. This simple and robust defi- 
nition is certainly meaningful for parabolic curves, but the probabilistic interpretation becomes 
ambiguous if there are secondary minima or other structures. 

An alternative possibility is to assume the validity of a certain CI model and to calculate 
the probability to observe the measured value for a given model parameter. Here the statistical 
error entering the x 2 function is taken from the model prediction. This approach is completely 
equivalent to using the log-likelihood function, and it has been verified that both methods lead 
to the same results. Each data point gets a varying weight depending on the model parameter, 
but independent of fluctuations in the data. For the present data this leads to more symmetric 
situations and an unbiased evaluation. In general, the resulting x 2 curves show a stronger sensi- 
tivity to details of the CI model, e.g. interference patterns, as can be seen in figure 5. The widths 
of the distributions are often, but not always, wider than in the previous case. The problem to 
extract limits from these distributions persists. 

The final limits presented here have been determined by applying a frequentist approach. 
Starting from a specific model with a scale parameter A true the cross section da/dQ 2 is cal- 
culated and then smeared according to the statistical error given by the predicted number of 
events. Distortions due to all uncorrelated and correlated systematic uncertainties (except for 
the parton distributions) are included, assuming Gaussian behaviour of the errors. The Monte 
Carlo experiment is then analysed in the same way as the data, i.e. with statistical errors from 
the prediction and including all systematics, resulting in a fitted value Af it . This procedure 
is repeated numerous times and the fit results are recorded in a probability distribution. The 
scale parameter is then varied and the 95% confidence lower limit A + (A - ) is defined as that 
value A true where 95% of the Monte Carlo experiments produce values of Af it which are larger 
(smaller) than the parameter A actually obtained from the data. Examples of CL distributions 
are shown in figure 5. For the present analysis the frequentist approach provides in general 
slightly weaker and more symmetric limits than the statistical method used previously [1]. 
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Table 1 : Lower limits (95% CL) on compositeness scale parameters A ± for various chiral struc- 
tures. Results are given for the present analysis of e~p and e + p data and for a combined analysis 
including e+p data at y/s = 301 GeV [2]. 
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Table 2: Coupling coefficients rf ab , fermion number F and 95% CL lower limits on M L q/\ for 
scalar (S) and vector (V) leptoquarks. Results are given for the present analysis of e~p and e+p 
data and for a combined analysis including e+p data at y/s = 301 GeV [2]. L and R denote the 
lepton chirality. The subscript 1 = 0, 1/2, 1 is the weak isospin. S and V differ by two units 
of hypercharge from S and V respectively. Quantum numbers and helicities refer to e~q and 
e~q states. 
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Table 3: Coefficients e q ab and 95% CL lower limits on Mq/X' for R p violating couplings to 
squarks. Results are given for the present analysis of e~p and e + p data and for a combined 
analysis including e + p data at y^s = 301 GeV [2]. 
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Table 4: Lower limits (95% CL) on the gravitational scale M s assuming positive (A = +1) and 
negative (A = — 1) coupling from the present analysis of e~p and e+p data and from a combined 
analysis including e+p data at y/s = 301 GeV [2]. 
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Table 5: Upper limits (95% CL) on the quark radius R q assuming point-like leptons (f e = 1) 
or common form factors (f e = f q ) for the present analysis of e~p and e + p data and from a 
combined analysis including e + p data at y/s = 301 GeV [2]. 
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Figure 1: Cross sections da/dQ 2 at y/s = 319 GeV for e~p — > e~X scattering (top) and 
e + p — > e + X scattering (bottom). HI data are compared with Standard Model expectations 
using the CTEQ5D parton distributions. The errors include only statistics and uncorrelated 
experimental systematics. Normalisation uncertainties are 1.8% (e~p data) and 1.5% (e + p data). 
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Figure 2: Fit results on the parameters e/A 2 of compositeness models using the combined 
e + p and e~p data. The inner and outer error bars represent one and two standard deviations, 
respectively, obtained by using experimental statistical and systematic errors. 
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Figure 3: Exclusion regions and lower limits (95% CL) on compositeness scale parameters 
for various chiral models obtained from the combined e + p and e~p data applying a frequentist 
method. 
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Compositeness: VV 
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Figure 4: NC cross section da/dQ 2 at \/s = 319 GeV normalised to the Standard Model ex- 
pectation. HI e~p and e + p scattering data are compared with curves corresponding to 95% CL 
exclusion limits obtained from each data set and from the combined data for VV composite- 
ness scales A + and A~ (top), couplings M/X of leptoquarks S£ and Vf (center), gravitational 
scales Ms assuming positive (A = +1) and negative (A = —1) couplings (bottom). The errors 
represent only statistics and uncorrected experimental systematics. 
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Figure 5: Examples of different methods to calculate limits on the compositeness scale A for 
VA, AA and LL + RR models obtained from combined fits including all data. Distributions 
°f X 2 ~ Xmin versus e/A 2 using statistical errors from the experiment (left) and statistical errors 
from the prediction (center). Confidence level versus e/A 2 from the frequentist method with 
dashed lines indicating the 95% CL limits for positive and negative interference (right). In all 
cases systematics are included. 



15 



